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quarks of generic mass value, we make use of a calculational scheme exploiting 
a large mass expansion together with a small mass expansion. The large mass 



. expansion for the effective action is given by a series involving higher-order 



Seeley-DeWitt coefficients, and we carry this expansion up to order l/{mp)^, 
where m denotes mass of the quark and p the instanton size parameter. For 
the small mass expansion, we use the known exact expression for the parti- 
cle propagation functions in an instanton background and evaluate explicitly 
the effective action to order (mp)^. A smooth interpolation of the results 
from both expansions suggests that the quark contribution to the instanton 
tunneling amplitude have a relatively simple m/o-dependent behavior. 
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I. INTRODUCTION 



Instantons as localized finite-action solutions of the Euclidean Yang-Mills field 

equations, describe vacuum tunneling and are believed to have important nonperturbative 
roles in low energy QCD. For an excellent review on instantons in QCD and general gauge 
theories, see Refs. [^Q. For actual instanton calculations, one needs to know above all the 
one-loop tunneling amplitude or the Euclidean one-loop effective action in the background 
field of a single (anti-) instanton. The latter quantity is thus of fundamental importance in 
instanton physics, and in the zero mass limit of scalar or quark fields 'tHooft |^ was able 
to calculate the appropriate one-loop contribution exactly. But, with finite quark mass, 
such exact calculation does not look feasible and one has to be satisfied with approximate 
results. [Note that, aside from up and down quarks, all other quarks possess sizable mass]. 
In this paper we shall describe our approach to determine the quark mass dependence in 
the one-loop vacuum tunneling amplitude, and report some new results from this analysis. 

By studying the field-theoretic effective action one can take systematically quantum na- 
ture of the fields into account, and already at the one- loop level it has provided us with 
certain relevant information on various physically significant effects 0]. In particular, the 
leading-order renormalization group coefficients in field theories are encoded in the diver- 
gences of the corresponding bare one-loop effective action. These divergent terms can be 
found most simply with the help of the background field method and the Schwinger- 
DeWitt proper-time algorithm [^,0]; they are entirely given by the second Seeley-DeWitt 
coefficient 02 [|7|,p|,p!0| in four- dimensional space-time. (For a recent literature discussing this 
method, see Ref. [llT|). But the evaluation of the full finite part of the one-loop effective ac- 
tion in any non-trivial background field corresponds to a formidable mathematical problem 
in general. For an approximate calculation (in a slowly-varying background) the so-called 



derivative expansion of the effective action has been utilized by various authors [12,13 



As for the contribution to the instanton one-loop effective action by spin-0 or spin-1/2 
matter fields of, say, mass m, we shall consider both the approximation applicable for rela- 
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tively large mp (p is the instanton size), i.e., large mass expansion and the mass perturbation 
scheme useful for relatively small mp. Note that the nature of the approximation is gov- 
erned by the dimensionless parameter mp. (Dependence on the renormalization mass scale 
p can be treated separately). The large mass expansion is essentially a series involving 
higher-order Seeley-DeWitt coefficients, for which a simple computer algorithm has been 
developed recently |T^-[l6|. We then make a smooth interpolation of the results found in 



those two different regimes, with the expectation that some general pattern, which is mean- 
ingful over a wide range of mass values, may emerge. This information should be valuable 
in phenomenological studies related to instanton effects. To connect the amplitude given 
for different mass scales, one should be careful about possible large finite-renormalization 
effects and renormalization schemes used. In this paper we treat various issues related to 
this general idea in a reasonably self-contained manner. 

In Sec. II we present a concise review on the Schwinger proper-time representation of 
the effective action, various renormalization schemes, and the large mass expansion. Also 
discussed are finite renormalization effects specific to renormalization prescriptions chosen, 
since they can introduce additional mass (as well as renormalization scale) dependences into 
the effective action. This understanding will become important when one has to change the 
results obtained in one renormalization prescription to that in another prescription. 

In Sec. Ill the one-loop effective action for a massive scalar field in a constant Yang-Mills 
field background is considered to see how our general scheme would fit in for this simple 
case. Here we make a detailed comparison between the known, exact, effective action (given 
in a single integral form) and the corresponding result based on the large mass expansion. 

In Sec. IV the spin-0 one-loop effective action in a Yang-Mills instanton background is 
studied on the basis of the large mass expansion. (Contributions due to fields of different 
spin can be related to this spin-0 amplitude). We consider up to the sixth Seeley-DeWitt 
coefficient. Here our finding is that, for mp ^ 1.8, the large mass expansion appears to give 
a good approximation to the effective action. 

In Sec.V we study the spin-0 instanton effective action for small mp, utilizing the known 



expressions for the massless propagators in an instanton background and the mass 



perturbation. Since the naive mass perturbation leads to a logarithmically divergent integral 
p!8| , a suitably modified perturbation method must be employed to obtain a well-defined 
small-mass correction term. We here reconfirm the (9((mp)^lnmp) term previously found 
in Ref. |T^, and provide for the first time the full 0{{mpY) contribution to the instanton 



effective action. 

In Sec. VI we consider an interpolation of our amplitude to intermediate values of mp, 
given the results of the previous two sections. Here we also make appropriate changes in our 
results so that they may describe the spin-1/2 instanton effective action; this result is directly 
relevant for quarks with nonzero mass. Note that, due to the hidden supersymmetry in an 
instanton background, one can utilize the result for the spin-0 case to find the contribution 
due to spin-1/2 fields 0. 

In Sec. VII we conclude with some remarks. In Appendix A some explicit expressions for 
higher-order Seeley-DeWitt coefficients can be found. Appendix B contains an analysis of a 
certain function which figures in our small-mass expansion of Sec.V. 

II. THE ONE-LOOP EFFECTIVE ACTION, RENORMALIZATION, AND THE 
LARGE MASS EXPANSION 

To be definite, we will consider a four-dimensional, Euclidean, Yang-Mills theory with 
matter described by complex scalar or Dirac spinor fields of mass m. Then, in any given 
Yang- Mills background fields A^^{x), one may represent the (Pauli-Villars regularized) one- 
loop effective action due to matter fields by 



T{A) = Aln 



Det(G-i + m2) Det(G'o ^ + A^^ 



(2.1) 



Det(Go ' + m2) Det(G-i + A2)_ 

Here, A = 1(— |) for scalar(spinor) fields, A is the large regulator mass, stands for the 
appropriate quadratic differential operator, viz.. 
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, —D^, (for scalar) 

G' = { ^ ' (2.2) 

(7/^)^, (for spinor), 

and Gq^ = G"^ 1^,^=0 = —d^. [Also, = D^D^ and = with the covariant 

derivative = — iA'^^^T"' = — iA^ ( denote the group generators in the matter 
representation satisfying the commutation relations [T", T^] = ifabcT'^ ), and our 7-matrices, 
which are antihermitian, satisfy the relations {'y^,'yu} = —26^^]- 
In the proper-time representation one represents T{A) by 

T{A) = -X r —(e'""'^ - e'^'')Ti[e-^'^'' - e^^^o (2.3) 
Jo s 

where 'Tr' denotes the trace over space-time coordinates and all other discrete indices. More 



explicitly, writing Tr = / d'^xtr and introducing the proper-time Green function (xs 



yj, r(y4) can be expressed as 
T{A) = -X r —{e-""'^ - e-^'') I (fxti (xs x) - (xs x)U=o , (2.4) 



The full effective action is thus determined if the coincidence limit (i.e., y = x) oi the 
proper-time Green function is known. The expression ( |2.4{ ) diverges logarithmically as we 
let A ^ CX3; to isolate such divergent pieces, we may exploit the asymptotic expansion 

0+ : (xs 



.n=0 



where the leading coefficient has the coincidence limit ao(a;, x) = 1. Using this expansion 



in (|2.4|) , we then see that the divergences in T[A) as A ^ 00 are related to the coincidence 
limits ai{x) = trai(x, x) and a2{x) = tra2(x,x), which correspond to the first and second 
Seeley-DeWitt coefficients respectively. Simple calculations yield 

ai(x) = 0, (for both scalar and spinor) (2.6) 

-j^tr(F^^(x)F^^(x)), (for scalar) 
a2(x) = < (2.7) 



tr(F^,,(x)F^i.(x)), (for spinor). 
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where F^y = F^JT^ = i[Dfj,,Dy]. (The 'tr' here refers to the trace over gauge group rep- 
resentation indices only). Based on these, we may now write the above effective action for 
7^ as 



T{A) = Y(ln 



A2 



with 



Y 



1 c 

12 (47r)2 ' 

1 c 



T{A) 



(2.8) 



3 (47r)2 



(for scalar) 
(for spinor), 



(2.9) 



(C is defined by tr (T"T ) = 5abC), and then the contribution 



T{A) = -X 



°° ds 



d X 



1 



s— + -s^ 
ds 2 ds"^ ' 



U=o 



tr 



X 



(2.10) 



becomes well-defined as long as is non-zero. [In ( p. 101) , (1 + s-^ + |s^^)|s=o/(s) 
/(O) + sf'iO) + \s^f" 



The logarithmic divergence in T{A) is canceled by the renormalization countert- 
erm associated with the coupling constant renormalization of the classical (bare) action 

J d'^xF^^F^j^. But the resulting renormalized one-loop amplitude depends on the renor- 
malization prescription chosen. From the very structure exhibited in (|2.8|) , our amplitude 
T{A) can be considered as defining a renormalized one- loop effective action; but, this pre- 
scription cannot be used for the strictly massless case. Instead, one may here consider adding 
to T{A) the counterterm 

(/i is an arbitrarily introduced renormalization mass) to obtain the renormalized one-loop 
effective action 

,2 



Ar(A) = -F(ln^) jd'xF;^F;, 



(2.11) 



TUA) = -Y{\ri—) d'xF^F- + T{A) 



12- 



(2.12) 



where T{A) is defined by ( p.lO| ). It should be remarked that ri.en(^), given by ( p.l2| ), is 
expected to have a well-defined limit for (i.e., does not exhibit infrared singularities), 

if the operator does not allow any normalizable zero eigenmode. 
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Other renormalization prescriptions may also be chosen. Let FmsI^) denote the renor- 
mahzed amphtude in the so-called minimal subtraction scheme [0] associated with the 
dimensional regularization. Then, to obtain the expression for FmsI^) scalar case, 

one should add to that of Frenl^) the following finite renormalization term 

C 



(2.13) 



(47r)2 ■ 12 

where 7 = 0.5772 ■ ■ ■ is the Euler's constant. In the thus found expression of FmsI^) ^^e 
mass parameter /i, which enters the expression through Fren(^), describes the normalization 
mass conventionally introduced in the minimal subtraction scheme. As for the expression 
of Fms(^) in the spinor case, the finite renormalization term to be added to that of Fi.en(^) 
turns out to be 



C 



;(ln4vr-7) j d'xF^^F^,. 



(2.14) 



(47r)2 . 3 

[This is the case when the spinor trace of 1 is taken to be four]. In another often- used 
prescription, one specifies the renormalization counterterm via the momentum-space sub- 
traction scheme, i.e., by imposing a normalization condition at certain external momentum 
value, p2 = ^2 'jj-^g^ corresponding renormalized expression, Fmom(^), is given by that 
of Fi.en(^) plus the finite renormalization term 



C 



(47r)2 • 12 



h In 



m 



3 



In 



'1 + 4^ + 1 



a + 4^ - 1 



or 



C 



(47r)2 . 3 



/ d^rF"- F°- 



5 , m m 
h In — + 4— 



(2.15) 
(for scalar) 



X / d'^xF^^F^^,, 




with the corresponding reinterpretation of the parameter /i. These renormalizat ion- 
prescription dependences of the one-loop effective action are of course explained by the 



fact that the tree-level contribution involves, as a multiplicative factor, the renormalized 
coupling \ (whose value may vary with renormalization prescriptions). 

Thanks to the exact connection formulas we have described above, knowledge on the 
one-loop effective action in one renormalization prescription can immediately be changed 
into that in another prescription. In fact, in theories containing several matter fields of 
different mass scales ( e.g., QCD with quarks of very different masses), one may well adopt 
different renormalization prescriptions for different matter field loops. We here note that 
use of the minimal subtraction for a heavy-quark loop is rather unnatural, due to the lack 



of manifest decoupling But this is not an issue in our discussions. 

The next task is to find the actual full expression for the one-loop effective action — 
at present, this is possible only with a background field of very special character. But, if 
the mass parameter is sufficiently large, it can be studied for generic smooth background 
fields by utilizing a systematic large-mass expansion, which is obtained by inserting the 
asymptotic expansion ( |2.5| ), say, into the formula ( p.lO[ ) for T{A). This assumes the form 



T{A) 



(m 



{an{x) = tra„(x,x)). 



[2.17) 



That is, for large enough mass, we have the one-loop effective action (in any renormalization 
prescription) expressed by a series involving higher-order Seeley-DeWitt coefficients dn{x) 



{n > 3), the calculation of which may be performed using a computer [^J-|T^. If only the 
leading term is kept with the series ( p.l7|) , one find, explicitly. 



m 



oo : T{A) = 



b^mJd'x tr[l{D,F,x){D^F,x)-^tF^uF,xFx^l 



16 

1 1 2 
327r2 m? 45 



(for scalar) 
(for spinor) 
(2.18) 



where D^F,^x = [D^,F,^x]- If the background fields under consideration satisfy the classical 
Yang-Mills field equations, one can show using the Bianchi identities and the property 
of trace that / d'^x tT[{D^F^x){D^F^x)] = 4:i J d^x tT[F^yFyxFxij\- Hence, for the on-shell 
effective action, (|2.18|) can be further simplified as 
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_L ]_ j_ 

n-2 ™2 



Jd'^x ii{F^yF^xFx^), (for scalar) 



m — > OO : T[A) = <^ (2.19) 

it^^iE / tr(F^^F^Ai^A/.), (for spinor). 

For some explicit expressions of the higher-order Seeley-DeWitt coefficients, see Sec. IV and 
Appendix A. Also note that the large mass expansion for the effective action in other renor- 
malization schemes can be obtained from the expansion ( |2.17] ) for T{A) and the exact con- 
nection formulas. 

The large-mass expansion is only an asymptotic series, and the useful range of the series 
( p.lTp (as regards the magnitude of m) will depend much on the nature of the background 
field and also on some characteristic scale(s) entering the background. For a sufficiently 
smooth background, this large-mass expansion may be used to obtain a reliable approxi- 
mation to the effective action even for moderately large values of m. But the series ( p.l7|) 
is bound to lose the predictive power for 'small' values of m, and for the small-m effective 
action one should employ a totally different strategy, such as the small-mass expansion if its 
exact expression in the massless limit has been known by some other methods. In the next 
section, we shall first see how good the large-mass expansion can be for the much studied 
case of the one-loop effective action in the constant Yang-Mills field background. Also con- 
sidered is its small-mass expansion which may serve, together with the result of the large 
mass expansion, as a basis to infer the behavior of the effective action for arbitrary mass. 



III. THE SPIN-0 EFFECTIVE ACTION IN A CONSTANT SELF-DUAL 
YANG-MILLS FIELD BACKGROUND 

In this section, various approximation schemes to be used later will be tested against the 
exact result, choosing a rather simple background field. In non-Abelian gauge theories, a 
constant field strength is realized either by an Abelian vector potential which varies linearly 
with x'^ or by a constant vector potential whose components do not commute In 



this paper we only consider the case of the Abelian vector potential. Assuming the SU(2) 
gauge group, an Abelian vector potential can be written as = —\f^uXuT'^ (with the field 
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strength tensor F^^ = f^,yT^/2), where is the third Pauh matrix. If we further restrict 
our attention to that with the self-dual field strength (i.e., F^^, = ^e^uxsFxs), its nonzero 
components may be specified by setting /12 = f^^^ = H with the constant 'magnetic' field 
H. 

In this Abelian constant self-dual field, let us consider the one-loop effective action due 
to an isospin-1/2, spin-0 (complex- valued) matter field, taking the mass m of our spin-0 
field to be relatively large so that the large mass expansion (|2.17| ) may be used. For this 
case, some leading Seeley-DeWitt coefficients are easily evaluated (using the formulas given 
in Appendix A, for instance). 



(3.1) 



4 



189 



{H/2Y 



as 



675 



{H/2Y 



Note that we get zero for all odd coefficients here. Using these values, we then find that, for 
relatively large m, the effective action is given by the series 



r{H;m) 



167r2 



1 H. 
120 



H 



H 



504 m^' 



720 m2' 



+ 



(3.2) 



where V denotes the four-dimensional Euclidean volume. 

For this case, it is actually not difficult to find the exact expression for the one-loop effec- 
tive action, following closely Schwinger's original analysis in QED After some algebras, 
one finds the trace of the proper-time Green function to be given by |p2|] 



tr {xs\x) 



(Atts) 



{Hs/2) 



sinh.\Hs/2) 



(3.3) 



One can easily check that the expressions given in (|3.1|) are correct ones by considering a 
small-s series of this exact expression. Inserting ( p.3|) into the formula ( p. 10 ) then yields the 
exact expression 



V(H;m) 



-2V r -e-^^-L 
JO s [Aits) 



smh\Hs/2) 3^ ' \ 



(3.4) 
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FIG. 1. Plot of the effective action T(H\m) 



Comparing the result of large mass expansion in ( |3.2| ) against this exact result, we can 
investigate the validity range of the former. From the plots in Fig.l, it should be evident 
that for mass values in the range m/\fH ^ 1, summing only a few leading terms in the 
series ( |3.2|) already produces the results which are very close to the exact one. The large 
mass expansion is useful if m ^ \fH . 

Now suppose that the exact expression ( |3.4D were not available to us. For mass value 
not larger than \/H, the large mass expansion ( p. 21) fails to give useful information. 
Nevertheless, if one happens to know the one-loop effective action for small mass, this 
additional information and the large mass expansion might be used to infer the behavior of 
the effective action for general, small or large, mass. [Note that, in an instanton 
background, this becomes a real issue since the full m-dependence of the effective action is 
not known there]. In exhibiting this, r{H;m) will not be convenient since it becomes 
ill-defined as m ^ 0. So, based on the relation ( |2.12| ), we may consider the renormalized 
action T^en{H;m, fi) given by 
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T,^n{H; m, jj) = - ) + T{H; m). 

(47r)^ • 6 fi'^ 



(3.5) 



which is well-behaved for small m. Large mass expansion for r^e^{H]m, fi) results once if 



the expansion (|3.2|) is substituted in the right hand side of (|3.5| ). 

To find the small-m expansion, we find it convenient to consider the quantity 

= T^cn{H; m, fi) - T^cniH; m = 0,fi) 



(3.6) 



In (|3;|), from (0) and 



167r2 



^4) 



-2C'(-1) 



(3.7) 



where C{s) is the first derivative of Riemann zeta function and ('{—1) ~ —0.165421. Notice 
that Q{H] m) is independent of the normalization mass fi and is well-behaved in the small 
mass limit. Explicitly, it is given by the expression 

{Hs/2r 



pm poo 

Q{H; m) = 2V drfi^ / ds 
Jo Jo 



(Atcs) 



smh\Hs/2) 



(3.8) 



and in the small mass limit, this leads to 

YH2 ^m^/H 



Q{H;m) 



87r2 
VIP 
167r2 



dm 



1/2 -m^ (log + 7) + ■ ■ 



-myn - {m^/H f{\og{m^/H) - 1/2 + 7) + ■ 



(3.9) 



In Fig. 2, graphs for Q{H; m), the exact one and those based on approximation schemes, 
are given as functions of X = m/\/H. The exact result, i.e., that based on the expression 
( p.8| ) is represented by a solid line, which exhibits a monotonically decreasing behavior 
starting from the maximum at X = 0. Clearly the small mass expansion up to 0{m'^ / H'^) 
provides a reliable approximation for X ^ 0.4, while the large mass expansion for Q{H] m), 



167r^ V 6 H 120 m"' 504 m"' 720 m"' 



(3.10) 
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(this formula is obtained from ( |3.2| ), ( p.5|) and (pl6|)) can be trusted in the range X > 1. 

In the intermediate region 0.4 ^ X ^ 1 the large- mass expansion curve (a long dashed line 

in Fig. 2) may then be smoothly connected to that given from the small-m expansion ( p. 91 ), 

assuming a monotonic behavior (as should be reasonable for a simple background field). 

Evidently, with this interpolation, one could have acquired a nice overall fit over the entire 

mass range even if the exact curve were not known. We also see from Fig. 2 the typical 

behaviors which are shown by the small-mass or large-mass expansion curves. 
Q/(V h' 

Exact result 

Small mass expansion up to 0(m7H) 
Small mass expansion up to 0(m^/H^; 
Large mass expansion up to Sg 



-2.0x10 




-4.0x10" 



0.5 1 1.5 

X (=m/H"') 



FIG. 2. Plot ofQ{H;m). 

IV. LARGE MASS EXPANSION FOR THE SPIN-0 INSTANTON EFFECTIVE 
ACTION 



We now turn to the case of a BPST instanton background i.e., a self-dual solution 
of Yang-Mills field equations given by 



(4.1) 



2 + 

where rifj_ya (« = 1,2,3) are the so-called 'tHooft symbols |Q and p denotes the size of the 
instanton. The associated field strength F^^, is 
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In this instanton background, the exact expression for the one-loop effective action due to a 
spin-0 or spin-1/2 matter field of nonzero mass is not known; only the result in the massless 
limit is known [^. This quantity will be studied with the help of approximation schemes in 
this paper. Specifically, taking the matter field to be that of an isospin-1/2, spin-0 particle, 
the corresponding effective action is studied using the large mass expansion in this section 
and by the small mass expansion in the next section. In Sec. VI, we then use these results for 
a spin-0 matter field to obtain the corresponding results appropriate to a spin-1/2 matter 
field (i.e., quark). Note that, in the case of a spin-1/2 matter field, a direct application of 
the small mass expansion can be very subtle due to the presence of normalizable zero modes 



for the massless Dirac equation p3|). 

The large mass expansion for the spin-0 effective action is described by our formula ( p. 17] ). 
To use this formula, one needs to know some higher-order coefficients in the series ( |2.5| ), with 
G^^ = —D^ and the instanton background given above. Calculations of these higher-order 
Seeley-DeWitt coefficients are straightforward in principle, but get very involved as the order 
increases. Fortunately, thanks to the rapidly growing computer capacity to handle a large 
number of terms in the symbolic calculations, the explicit expressions for the Seeley-DeWitt 
coefficients in general background fields have been found recently up to the sixth order 



14-16 . We will utilize these results for our calculations below. 



In the instanton background ( [4.1|) the renormalized one-loop effective action f (A), de- 
fined by ( |2.10| ), will be a function of mp only. Hence our large mass expansion is really an 
expansion in 1/m^p^. Also the expressions for the Seeley-DeWitt coefficients are simplified 
considerably if we take into account the fact that our background field satisfies the classical 
Yang- Mills equations of motion. For such on-shell background fields, the space-time integral 
of the Seeley-DeWitt coefficients a„(x),n = 3,4,5 (for a spin-0 matter field) are given as 

n 

J d^xdsix) = ^Jd^x tr [F^^F.pFpf,] , (4.3) 
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2 

35' 



17 

2 -j^g FpuFfjLvPxuFxK ~\~ '^^FpvFupFpxFxp 



420 ^f^'^'^p'^'^^^'^'^P°' 



105 ^'^^^^P^p'^^'^^^ 



(4.4) 



y (i^x a5(x) = y d^x tr 



7 p F F F F 

945 '^'^ ''"'^ '^^ 



47 

p F F F F 

^ j^2g '^'^ '^'^ ''''' 



F F F F F -^i—F F F F F F F F F F 

izD Do ioy 

37 4 2 

'^^'^^^F^iyFupFpuFo-YF-j-^ -\- Fi/t-Fyct i^F) pFyp) (^F) ^Fpu^ gg F\i^(^F) p^Fup^Fi,p(^F) ^Fx,^ 

2 4 2 

— -r^Fxu{Dp_Fyp)Fpfj{DpF^\) + —F„rFur{F> pFyp){D p_Fj,p) + —Fp_rFrcT{F>pFyp){DfjF^p) 
ioy 00 00 

4 „ „ „ J ^^^^^ 



189 



F(jtFyi'(^F)^Fup) {D pF^ 



per) 



Note that the on-shell expressions for the space-time integral of d^ix) and 04 (x) involve only 
the field strength, while that for 05 (x) involves the derivatives of the field strength also. For 
the expression of di^ix), which occupies more than a page, see Ref. [1^. In Appendix A, the 
expressions valid without using the classical equations of motion (and before the space-time 
integration) can also be found. 

Inserting the expression ( [4.2| ) for the field strength into the formulas ( ^.3| ) and 
carrying out tensor algebra and trace calculations, we find 

64p6 16 



and 



y d^x d-i(x) = J d 
d'^x 0,4 (x) = / d'^x 



15(x2+p2)6 75 p2' 

,4^ 544p8 _ 272 TT^ 
^ 35(x2 + p2)8 - 735^- 

The next coefficient 05 involves the covariant derivative of field strength. 



(4.6) 
(4.7) 



D.F 



X^ pu 



(4.8) 



(x2 + p2)3 

Calculations of higher-order Seeley-DeWitt coefficients with the instanton background can 
be very laborious. Together with the formulas given above and that in Ref. ||T5[ for a6(x), 
we have thus used the "Mathematica" program to do the necessary trace calculations as 
well as tensor algebra. From the expression for the 05 coefficient, we obtain the result 
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J d^x {x) = J d^x 



512(35a;V - 39pi°) 



2835 p6 ' 
while, for the dQ{x) term, 

256 



315(x2 + p2)io 
1856 

(4.9) 



j d^x dG{x) — — J d 



51975(p2 + a;2)i2 

X [397710xV^ - 765270a;V^ + 404961a;V - 86418a;V^° + 2876^^^] 

= ^^4. (4.10) 

444675 p8 ^ ^ 

Based on the exphcit calculations given above, we obtain the following large-mass ex- 
pansion for r(mp): 

- 1 1 17 1 232 1 7916 1 

[mp) - -T^^^ - + ~ 148225 m^p^ ^ ' ^ 



r(mp) 



-0.01 




Large mass expansion up to 
Large mass expansion up to a. 
Large mass expansion up to §5 
Large mass expansion up to a. 



mp 



FIG. 3. Plot of T{mp) for the instanton background. 



In Fig. 3 we have given the plots based on this expression (first keeping only the fis-term, 
then including the a4-term also, etc). This is a useful approximation when mp is large, say, 
mp ^ K. What would be the lower-end value K here ? In the absence of the exact 
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expression for T{mp), a possible criterion for telling the validity range of the series ( [4.11|) 
will be as follows. If Ai denotes the 0{ ) term in the series and Vi = J2n=i An, we may 
demand that the series ( [4.11[ ) remain stable in the sense that the relative importance of 



each newly added term decreases, i.e., 1 



> 



> 



> 



> ■ • 



As this criterion 



is used, we obtain the (conservative) value K ^ 1.8. The result of large mass expansion 
may thus be trusted in the range given by mp ^1.8. 



V. SMALL MASS EXPANSION FOR THE SPIN-0 INSTANTON EFFECTIVE 
ACTION 

For small mp, that is, mp significantly below 1, the one-loop effective action in the 



instanton background ( [4.1|) can be studied with the help of the small mass expansion or mass 
perturbation, since its exact expression in the massless limit is known. Here we shall denote 
the corresponding spin-0 effective action Fren, which is defined by ( p.l2| ), as Fren(^, P, /^)- 
For m = we have, from the computations of 'tHooft 0, 

Fren("^ = 0, p, yu) = ^ In /ip + a(l/2) (5.1) 


with a(l/2) = I7 + |ln7r — ^C'(2) ~ fi — 0.145873. Our goal in this section is to com- 
pute explicitly the (9(m^) term of Ty.f,^{m, p, p). Note that this small-m approximation for 
^ren{^, P, fJ') coutaius a non-analytic piece in m and so it is not a trivial task to extract the 
desired term. 

Our first task is to develop a small mass expansion for Fren(^, P, p), which is finite at 
every order. For the purpose it is convenient to consider its derivative with respective to 



m 



2 



i.e., dTj-en/dm'^ , which is independent of the normalization mass p. The latter, being 



equal to the m^-derivative of our regularized effective action T[A), will have the proper-time 
representation 

dT 

ren 

(m, p, p) 



dm"^ 



lim 



p ds e-™'^ 


j d'^x tr 




— (^xs 


x) 










1 A^=oJ 



(5.2) 
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where we have used (12.41). Here note that Ixs 



x) = hniy^x (xs 



y) is nonsingular as long as 



s > 0. Paying due attention to various (singular) limits involved, it is possible to recast the 
formula ( p. 21) into that involving ordinary spin-0 propagators 



G{x, y; m^) = (x 



y), Go{x,y;m^) = (^x 



— + 

The explicit formula, which can be derived from ( p. 2] ), reads 
dT 

rcn 



(5.3) 



dm? 



d X lim tr[G(x, y;m ) — ao{x, y)Go{x, y;m )], 

y 



(5.4) 



where ao(x,?/) is the zeroth order coefficient in the asymptotic series ( |2.5|) . For small {x — 
y)^ii O'oi.^iy) has the following expression 

ao(x, y)= I + i{x- y)^A^{y) + ^{x - y)^{x - y),,[df,A^{y) + d^A^{y) + i{A^{y), A^{y)}] 



+Oiix~y)'). 



(5.5) 



Presence of the a^GQ term in ( [5.4| ) guarantees a finite result for dTj.f.^{m, p, fi) /dm"^ . 

For small m, one may then try to evaluate the right-hand-side of ( p.4| ) by exploiting the 
appropriate expansion of the propagators in and the known exact massless propagation 



function in the instanton background. We shall denote the latter by G{x, y) = {x 
But a naive expansion of the form 



-D2 



G{x, y; m^) = G{x, y) — j d'^zG{x, z)G{z, y) + j d^zd'^wG^x, z)G{z, w)G{w, y) 

+ ••• (5.6) 

is not valid since, aside from the leading term, all other terms of this series involve loga- 
rithmically divergent integrals. (Note that G{x,z) = O(^) as |z| ^ oo). Moreover, the 
= limit of dV^cn/dm? does not exist in the instanton background since, according to 
explicit calculations (see the comment after ( |5.22| ) below), the integral in the right hand side 
of ( |5.4| ) for m? = diverges logarithmically. This indicates that, as 0, T^-eni^, p, p) 

approaches the 'tHooft result ( p.l| ) in a non-analytic manner. To resolve this problem, we 
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shall below describe an alternative expansion scheme (which utilizes the idea of Carlitz and 



Creamer |^ in a suitable form). 

The expansion we shall use has the form 



1 



1 



1 



E 

Lr=0 



-m 



(-D2) (-92) 



1 



(5.7) 



This can be derived in the following way. First observe that 

1 1 , 1 f 1 1 



—D"^ + m? — 52 + m? 
Then, using the identity 



1 



(5.8) 



1 



1 — m 



1 



it is not difficult to see that (15.81) can be rewritten as 



(5.9) 



-D^ + w? i-D 



— m 



1 1 



92) J —92 + m2 

(5.10) 



This last equation may be solved for 



-7^2^ — 2 in an iterative manner, and the result is the 



expansion ( ^.7|) . Evidently, ( |5.7| ) is an expansion in powers of m'^{ ~ (_ga) }(— ^^) ^ga+m^ ; 
and we expect that this will yield a convergent series for small m if the background field is 
such that 



X 



1 



1 



?/) — > (sufficiently fast), as \x — y\ 00. (5.11) 



^ -D2) (-52) ^ 

In the case of the instanton background, this means that we have to work with the expression 
given in the singular gauge, i.e.. 



x^(x^ + p"^j 



(5.12) 



[Here, rj^^a differs from rj^^^a only by the sign in the components with /i or z/ equal to 4] This 
is allowed by the gauge invariance of the effective action. 
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If only the leading term of ( |5.7| ) is used in ( |5.4| ), we now find that 



dV 

ren 
dTH? 



d X lim tr 

ao{x,y)(x 



X 



,2 



+ C(m"). 



(5.13) 



—(9^ + 

This is not yet in the convenient form for actual computations. So, based on the following 
observations 



X 



1 



y = {X 



y) + (x 



1 



1 



1 



-9^ + 



ao(x,y)(x 



— (9^ + 



;-a2) 



y) + {x 



+ (terms vanishing as y ^ x) 



y,. 

(5.14) 
(5.15) 



we make suitable rearrangements in the right hand side of ( ^.13| ) to obtain the following 
formula (to be used for computations) 



1 d^x 


lim tr < Ix 


1 









dTren{m,p,lj) 

dm? 



where the function J{x) is given by 

1 



y) -^o(a;,2/)(x|y-^|i/) j> - J{x) 



(5.16) 



J{x) 



d z ti < (x 



-92) 



X 



(5.17) 



But for the J{x) term, what we have in the right hand side of ( ^.16|) is just the (logarithmi- 
cally divergent) expression representing ^^^^|^^T^^|m2=o (see ( p.4|) ). As we shall see below, 
this divergence is tamed by the additional term J{x). The very structure of J{x) given in 
( p.lTp also ensures that it is free of any short-distance divergence. 

The first term inside the integrand of ( |5.16| ) is evaluated as follows. The spin-0 (and 
isospin-1/2) massless propagator in the instanton background ( |5.12D is given by |T7 



G{x,y) = (x 



1 



1 + 



p'^(x-y+ir]^ 



47r2(x — y)2 
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(5.18) 



Then, writing x = ?/ + e, we find after some straightforward calculations 



X 



1 + 



p {y ■ ^) - p y {y + p y 

2(^2)2(^2 + ^2)2 



\ y'^iy'^ + p'^) J y'^iy'^ + p'^) 



0{e) 



(5.19) 



On the other hand, if (|5.12|) is inserted into the expression (5^), we have 



ao{x,y) = I + i{x-y) 
1 



{x -y)^{x -y) 



p'^VfMuaTay„ 

"y'iy' + p') 

p^iy'^S^u - y^,yu) + 2ip^{2y'^ + p'^)f]^xaTayuyx 



(^2)2(^2 _^p: 



,2^2 



+ Oiix-y)') (5.20) 



and therefore 



tr <^ ao{x,y)(x 



p y — p {y ■ e) 

2(^2)2(^2 + ^2) 



+ 0{e) 



(5.21) 



From ( |5.19| ) and (|5.21|) , we thus obtain the following expression 



lim tr < ( X 



y) - ao{x,y)(^x 



A7T^{X^ + p2)2' 



(5.22) 



[We here remark that the result (|5.22|) is unchanged even if one takes the regular-gauge 
instanton solution (|4.1| ) as the background field]. Clearly, with this term alone, the remaining 
x-integration would yield a logarithmically divergent result. 
We now turn to the evaluation of J{x). Noting that 



( ^.17] ) may be rewritten as 



-92) 



X ) = — m 



-ip-{z—x) 



(27r)>2(p2 + ^2 



(5.23) 



J(x) = 



(27r)>2(p2 +^2) 



F{x,p) 



(5.24) 



with F[x,p) given by 

F{x,p) = J d'^ze-'P-'tTi^(x\ (^j^ 



2 / d^z e-'P-^-d, 



1 



p^X-Z 



^Tl\x-Zf \^^l+p2/3;2^1 + p2/^2 



(5.25) 
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In ( p.25|) we have used the expression (|5.18| ) and the factor 2 at front arose from the isospin 
trace. We are here interested in 0{1) or (9(logm^) contribution to the right hand side of 
( p.l6| ). Let us see when and where such contribution can arise, based on our formulas ( |5.24| ) 
and ( |5.25| ). For any finite x- value, the function F{x,p) is well-behaved for all p. Due to the 
overall multiplicative factor in ( [5. 241 ), J{x) for finite x (or, more precisely, for x'^ satisfying 
the condition |x| <^ ^) would then be 0{rn?) and hence no desired contribution. It is thus 
sufficient to study J{x) for large x, i.e., x^ in the region |x| > L with p <^ L <^ irT^. Now, 
due to the factor e^'''^ /p'^{p^ + m?) within the integrand of (|5.24|) , we further conclude that 
the small-p region of F{x,p), with |x| > L, can be the source for the desired contribution; 
if the contribution from the region \p\ ^ m is excluded from the right hand side of ( ^.241) , 
J{x) becomes 0{m?). 

To study the function F{x,p) for |x| > L (with p L m^^) and \p\ ^ m, we write 
F{x,p) as the sum of its value at p = plus the correction term, viz.. 



F{x,p) = F{x,p = 0) + C{x,p) 



(5.26) 



Then, from ( |5.25D , 

F{x,p = 0) = 2jd^z (-a 

1 



1 _|_ P a^-z 



lim 



..... , ^ d'QR'^ 

R^oo 271^ J\z\=R R 



+ 



(5.27a) 



1 d ( 1 + 4^ 

1 x^z^ 



(5.27b) 



where we used Gauss's law. [Note that, for very large \z\^ the integrand in ( [5.27aD behaves 
like 0(^)]. Evaluating the surface integral in (|5.27b| ) immediately gives 

1 



F(x,p = 0) = 2 



- 1 



and hence 



F(x,p = 0) = -^ + O(^) 



for Ixl > L. 



x^ 



x^ 



(5.28) 



(5.29) 
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Now, inserting the thus evaluated F{x,p = 0) for F{x,p) in (|5.24|) , we obtain the following 
contribution to J{x): 



J{x) 



1^1 -^'/"^l^-j;^)'"' (4) 

' [l-m\x\K^{m\x\)]e{\x\-L) + 0{m^). (5.30) 



(Note that we have assumed niL <^ 1). On the other hand, it is possible to show (see 
Appendix B) that C{x,p) in ( ^1261) is at most 0{\p\L^) or C(^) or O(f^), when |x| > L 



and \p\ ^ m. With this finding used in ( |5.24| ), it is easy to see that no 0{1) or 0{logm) 



contribution results from the C{x,p) part of F{x,p). Thus, to the order we want, our 
formula ( p.3(J| ) has no further correction. 

Evidently, if the contribution in ( |5.30| ) is considered together with that in ( |5.22| ), the 



x-integration in ([5. 161) will give a finite result. Furthermore, since the function F{x,p) does 



not involve mass m at all, the scale L we introduced can be chosen, for mp 0, such that 
p <^ L <^ 1/m. With this understanding, we may now perform the integral in the right 
hand side of ( ^.16|) to secure the unambiguous result 



^— ^ = — In(mp) + y(7 + - - ln2) + C(m p ). (5.31) 



Then, based on this formula and the 'tHooft result ( |5.1| ), we immediately obtain the desired 
small-mass expansion for Tj-e^^m, p, fi): 



T^en{m, p, p) = T^en{m = 0, p, /i) + / 

Jo 



^^2 9^rcn{m,p,p) 



= iln(pp) + a(l/2) + [In(mp) + 7 - In 2] + 0{{mpf). (5.32) 



The (9((mp)^ In(mp)) term in this formula was first found in Ref. [TH|, while the 0{{mpY) 
term without the In(mp) factor is new. 
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VI. MASS INTERPOLATION AND THE SPIN-1/2 INSTANTON EFFECTIVE 
ACTION 



In the previous two sections the spin-0 instanton effective action were computed for 
relatively large mp and for small mp. The result can be summarized by 

1 



Trenl"^, p, p) = - In(yup) + a(l/2) + Q{mp) 
o 

with the quantity Q, a function of mp{= X) only, behaving as 
iX2lnX + i(7-ln2)X2 + ..., 



_\__\_ IL J_ _|_ 232 1 _ 7916 1 

75 X2 735 X4 2835 148225 



+ 



(6.1) 

(X < 0.5) 
(X > 1.8). 

(6.2) 



(Note that, in the instanton background, ( |2.12|) implies that Vj-^^im, p, p) = ^\n{p/m) + 
r(mp)). In the indicated validity ranges of X, the function Q{X) is plotted in Fig. 4. We have 
here assumed that our small mass expansion in ( [5.32[ ) can be used reliably for X ^ 0.5; this 
estimate is based on measuring the effect of the {mpY term (with the numerical coefficient 
taken to be 0{1)) against the terms which appear explicitly in ( p. 321 ). 



Q(X) 



-0.1 



-0.2 



-0.3 





Small mass expansion 




Large mass expansion up to a^ 


\V 

\ \^ 


A plausible form for Q(x) 























0.5 1 1.5 

X(=mp) 

FIG. 4. Plot oiQ{mp). 
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Now, what could be said on the behavior of the function Q{X) in the intermediate 
region 0.5 ^ X ^ 1.8? Since the background field under consideration has a smooth profile, 
one naturally expects that Q{X) also be a smooth function of X; that is, Q{X) would be 
represented by a smooth interpolating curve connecting the known forms of the curve in 
the regions X ^ 1.8 and X ^ 0.5. Let us further assume that the region for interpolation, 
0.5 ^ X ^ 1.8, can be viewed as being reasonably small. Then, looking at how Q{X) 
actually behaves for X > 1.8 and X ^ 0.5 (see Fig.4), it appears to be quite plausible to 
suppose that Q{X) is a monotonically decreasing ivMciion of X for all X > 0. But, since we 
have in no way proved this monotonic behavior in the presence of the instanton background, 
one may regard this as a conjecture. [For instance, the possibility that Q{X) may develop a 
local maximum or minimum within the range 0.5 ^ X ^ 1.8 is not excluded. Incidentally, 
such monotonic behavior was also observed in the case of a self-dual constant field strength 
(see Fig. 2)]. Accepting the conjecture, it might be useful (especially for phenomenological 
analysis of instanton effects) to have a certain smooth function Q{X) in the entire range 
X > which meets this requirement. With the plausible curve for Q{X) taken by that 
given in Fig.4, we have found (after some trial and errors) that it may be described by the 
function of the form 

, , 1 ilnX + a-(3 a + /?)X2-ix4 ^ ^ . ^ 

with a = a;(l/2) ~ 0.145873 and (3 = |(ln2 — 7) ~ 0.05797. This form incorporates 
correctly the small-X and leading large-X behaviors shown in ( |6.2|) . For the tunneling 
amplitude which is more directly related to e"'""™, this amounts to using the expression 
(PD with 

g-r,en(m,p,M) = (;^p)-i/6g-a(i)-Q{mp) ^f^^ arbitrary mass m). (6.4) 

Various results obtained for the spin-0 field case can be used to derive the corresponding 
results appropriate to the spin-1/2 one-loop instanton effective action. The latter will be 
needed if one wishes to consider the loop correction to the vacuum tunneling amplitude in 
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QCD due to quark fields. In a self-dual Yang-Mills background, the hidden supersymmetry 
of the system allows one to express the spin- 1/2 proper-time Green function (xsly)^^/^-* = 
{x\e~'^^^^^'^\y) in terms of the corresponding spin-0 function {xs\y)^^^ = {x\e~^^~^'^^\y) (with 
the same isospin representation assumed). Explicitly, this is described by the operator 
relation 

e-^(-^)^ = e-^(-^^)i±^ + lD^.-<-^'hD^-^ + P, (6.5) 

z — iJ^ z 

where P is the projection operator into the zero mode subspace of and can be expressed 

by P = (1 - 7D-^7D)(i^) Using the relation (U) with the definition of V,^^{A) 

(see Sec. II), it is then possible to derive a simple relationship between the spin-1/2 and 

spin-0 one- loop effective actions. If r(^/2)(^) (r(°^(A)) denotes the one-loop effective action 

as defined by (|2.12|) for a spin-1/2 (complex spin-0) field of mass m in a self-dual Yang-Mills 



background, we have in fact 

1 



= --n^ln(^) - 2rW(A) (6.6) 



or, for the respective contributions to the tunneling amplitude 



py«,g2r-(A)^ (6.7) 



where is the number of normalizable spinor zero modes in the given background ^3 
Now, using the result ( |6.1| ) for T[^l^(A), we have the spin-1/2 instanton effective action 
expressed as (with np = 1) 

rili^\m,p,fi) = - In - - i In /ip - 2«(l/2) - 2Q{mp), (6.8) 

or, for the tunneling amplitude. 



where Q{mp) is the function specified in ( |6.2[ ) (and represented in Fig.4). 

The expression in ( p. 8] ) or ( |6.9| ) describes the one-loop contribution to the vacuum tun- 
neling by an isospin-1/2 quark field of mass m. If one accepts our conjecture, the function 
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Q{mp) may be taken as a monotonically decreasing function of mp which has the hmiting 
behaviors as given in ( |6.2| ). The renormahzation prescription appropriate to the expression 
( |6.8[ ) is that specified by ( p.l2| ). If one wishes to obtain the corresponding amphtude in the 
minimal subtraction in the dimensional regularization scheme, the finite renormahzation 
counterterm (see ( |2.14| )) 

^^-l-^(ln47r - 7) I d'xF^^^F^, = ^{InAn - 7) (6.10) 

must be added to the expression (|6.8| ). Thus, in the minimal subtraction scheme, the 
amplitude due to a spin- 1/2 quark of mass m reads 

/ \ 1/3 

:'iVJUm,p,u) ^ I \ ,2a(l/2)+2Q(mp) 



fi V47re-T/ ■ ^ ■ ^ 

With Qijnp) set to zero, this reduces to the result of 'tHooft 0. For applications to the real 
QCD with the SU(3) gauge group, one must also take into account the well-known group 
theoretical factor associated with various ways of embedding the SU(2) instanton solution 



VII. DISCUSSIONS 

In this work we studied the massive quark contribution to the one-loop instanton effec- 
tive action in QCD. For this purpose, we made use of the approximation scheme valid for 
relatively large mass as well as the small-mass expansion. These considerations provide a 
reliable approximation to the one-loop effective action if the magnitude of mp is such that 
mp ^ 1.8 or mp ^ 0.5. The expression for the effective action contains a function Q{mp), 
the magnitude of which is uncertain in the range 0.5 ^ mp ^ 1.8. Based on the known 
behaviors of Q{mp) in the ranges mp < 1.8 and mp ^ 0.5, we suggested that Q{mp) be 
a smooth, monotonically decreasing function of mp. If the latter turns out to be true, a 
simple interpolation formula for Q{mp) (as we considered in ( |6.3D ) suffice for considerations 
in most phenomenological analyses. 
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What can be done to reduce the uncertainty in the function Q{mp) for 0.5 ^ mp < 1.8 
? With the exphcit calculation of the 0{{mpY) term in the small mass expansion, it should 
be possible to push the lower end of the uncertain range to a slightly higher value. On 
the other hand, we expect that including the next higher Seeley-DeWitt coefficient in the 
large -mass expansion would not bring a significant new information. More useful direction 
might be to try a direct numerical evaluation of the functional determinant, with the help of 
the scattering theory in a radially symmetric background field. (Some related techniques are 
discussed in Ref. |^6|) Perhaps, by some mathematical argument, it might also be possible 
to actually prove that the function Q{mp), which is equal to ri.en("^, P, p) — Trenl^^ = 0, p, p) 
( for a spin-0 field) in an instanton background, is a monotonically decreasing function of 
mp. These are left for further study. 
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APPENDIX A: 



In the literature [1^-16], the Seeley-DeWitt coefficients a„(x) for a spin-0 or spin-1/2 
matter field have been calculated up to n = 6. Here, for the case of a (complex) spin-0 field, 
we shall give the explicit expressions for a„(x) up to n = 5 in a general off-shell background 
field. They read 



6 



tr 



.2 1 

TTT^KX^X^PfjivFuK ~l~ . (-,„ FnxPpivPxuFviJ, ~l~ „ r ^i^X^XK^piuPu^ 



(Al) 



21 



420 



35 



+^F^xFx^.F^,F,^ + i^F,x{D^Fxu){D^F,,) + i^F^x{DxF^,){D,F,^) 
+ ^{D^DxF^,){DxD^F,^] 



(A2) 



1 



tr 



945 



Ff^xFx^F^jjFypFpf^ i ^„Fi;\F\iiFt;,jFiinF, 



63 



nX^ X^^ Ku^ pu 



120 

.16 .22 .31 

^ , ^ FnxFpvFxpFpKFpp ~\~ ^ -| Q^Fi^xFxpFi^pFppFpp^ -\- i „ _„ Fi^xFxp,Fi^pFpi^Fpi, 



189 



378 



945 

53 1 1 

+i-z;i^Fi^xFxKFfj,uFupFpp + -FfixFxp{DuF^p){DyFpp) + i—Ft,x{F>xDpF^p){Dp_D^Fpi^j 
o fo 9 io 

--Lf^^Fx^(D^F,p){D^Fp,) + ^{D^DxDpF,p){DpDxD^Fp,) + ^F^^Fp,F,x{DxD^F, 
ioy ZdZ oIo 

+i^{D^Fx^,){D^D,F^p){D,Fpx) + ^F,x{DpF,p)Fx^{D^Fp,) + i^F^x{D^D,Fxp){D,D^Fp^) 



up) 



63 



63 



5 5 

-^Fkx{F> pFxu)Fkp{D pFpjy) + —Fi^x{F>p_F^p)Fpx{Dp_Fi, 



21 
5 

5 10 
+^Y^(^aFp.)(/^kFam)(/^m^kF.p) - —F^x{DxFp,)F,p{DpFp^) - ^Fp^F,A(I^AFp.)(/^.F«p; 

+^Y^(/^KFAM)(/^MF.p)pA/^.i^p.) - ^F^xiDpFxu)F^p{D^Fp,) + ^F^x{DpFx.)F,p{D^Fp,) 

+-rz7:Fi^x{F> pFyp)Fpy{D p^Fxk) — -z;i^Fp^F^xFxpiF>pDi^F^p) + —F,xFx.{DpF,p){D^ 

Fpv) 

ioy OIO ZoZ 

16 16 19 

— —Fkx{F> xFnu)F^p{D ^Fpn) — -——F^xFxnFup{DnDf^Fpjj) — ——Ffix{F>xFnu)F^p{DpFjjn) 
DO 189 756 

19 25 26 

— -;^Ff^x{F>pFj,p)Fpx{F>K,Fpy) + -^Fp^jFt,x{F> pFxu){D pF^p) — -^FpyFf^x{F> pFxu){D pFp^,) 

34 41 61 

-Y^^pf^^^^>-(^>-^t^^)(^ p^^^^) ~ '^^Ff^xFxp{DpF^p){Df^Fpy) + j^FpjjFf^x{D pFxk){D pFyp) 

« 1 ,^3, 



756 



F^xFpp {DpFx^){DpF,p) 



where Da^^p. = [Dx, Fp,] and D^DxFp, = [D^, [Dx, Fp,]], etc.. 
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APPENDIX B: 



The function C{x,p) in ( |5.26 ) is given by 



C{x,p) = 2jd'z{e-'P-^-l) i^-dl 



1 + ^ 



, (Bl) 



and we are here interested in its behavior for |x| > L(with p <^ L m^^) and p ^ m. We 
divide this quantity into two parts, i.e., C{x,p) = C^{x,p)+C^{x,p), where C^{x,p) denotes 
the contribution with the region of integration restricted to \z\ < Li(with p <^ Li < L) and 
C>(x,p) that from the region \z\ > Li.(We take Li to be of the same order as L.) Then, 
for C^{x,p), it will be safe to make an approximation e~*^'^ — 1 ^ —ip ■ z (i.e., 0{\p\L) at 
most) inside the integrand of ( [Bl| ) and so we find immediately 

^2 



C^{x,p) = 0{\p\L!-), (for|x|>L). 



x^ 



On the other hand, for C>(x,p), we may expand the factor 
of ( pT|) as a power series in p/\x\ and p/\z\, and then 



(B2) 

in the integrand 



1 / 1 + 



47r2(x — z^ 

.2 



y 1 + pV^y 1 + pv^2 



p 



3(x2)2 + 3(^2)2 ^ 2x2^2 _ 4^ ■ z{x'' + z'') 



in'^x'^z'^ 



47r2(x — 2;) 2 

¥2 



+ 



(B3) 



8(x2)2(22)2 

As the differential operator —dz acts on this expression, the first term in the right hand 
side of ( |B3|) can be dropped. Thus, we may write 



C>ix,p) 



2 dhU 

J\z\>Lx ^ 



-a 



-%p-z 



P^ 



1 



3(x2)2 + 3(^2)2 ^ 2x2^2 _ 4^ . ^/^2 ^ 



+ ■■ 



(B4) 



_47r2(x-z)2 8(x2)2(z2)2 

and, for |x| > L and \p\ ^ m, it is not difficult to show that this can only lead to terms of 
0{\p\Lj^) or O(^) or O(f^). Hence, C(x,p) is at most 0{\p\L^) or 0{^) or 0{^). 
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